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ferences therein). For planets with a significant magnetic field, duce reconnection leading to auroral brightenings (Cowley et al. 
namely Jupiter and the Earth, the X-ray emission has been found 2005 ); compared t0 Earthi though? where aurora brightenings are 
to separate essentially into two components: a high-latitude au- yery rapid ( „ tens of minutesX those on Saturn can last for days 
roral one, and a low-latitude disk component, which has different (clarke et aLE 005). Moreover, UV auroral activity has also been 
spectral properties, relating it to solar X-rays (Maurellis etal. obse rved on Saturn during a period (Oct. - Nov. 2005) of quiet 
2000 Gladston e et al |2002 Branduardi-Raymont et al. 12003 so i ar wind conditions, suggesting an intrinsically dynamic mag- 
Bhardwaj et al. |2005cj Eisner et_aL|2Q05l Branduardi-Raymont netosphere with in j e ctions of plasma occasionally taking place 
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is not the case. While the characteristics of Jupiter's UV aurorae they find evidence for a direct relationship between Saturn's au- 
are dictated by the rotation of plasma internal to its magneto- roral actiyity and solar wind conditionSi while the corre i at ion is 
sphere (Bhardwaj and Gladstone l2000l Clarke et al. 120041 1. those nQt SQ strong at j up jt e r 



Received accepted ... 

ABSTRACT 

Aims. We approach the study of Saturn and its environment in a novel way using X-ray data, by making a systematic and uniform 
spectral analysis of all the X-ray observations of the planet to date. 

Methods. We present the results of the two most recent (2005) XMM-Newton observations of Saturn together with the re-analysis 
of an earlier (2002) observation from the XMM-Newton archive and of three Chandra observations in 2003 and 2004. While the 
XMM-Newton telescope resolution does not enable us to resolve spatially the contributions of the planet's disk and rings to the X-ray 
flux, we can estimate their strengths and their evolution over the years from spectral analysis, and compare them with those observed 
with Chandra. 

Results. The spectrum of the X-ray emission is well fitted by an optically thin coronal model with an average temperature of 0.5 keV. 
The addition of a fluorescent oxygen emission line at ~0.53 keV improves the fits significantly. In accordance with earlier reports, 
we interpret the coronal component as emission from the planetary disk, produced by the scattering of solar X-rays in Saturn's upper 
atmosphere, and the line as originating from the Saturnian rings. The strength of the disk X-ray emission is seen to decrease over 
the period 2002 - 2005, following the decay of solar activity towards the current minimum in the solar cycle. By comparing the 
relative fluxes of the disk X-ray emission and the oxygen line, we suggest that the line strength does not vary over the years in the 
same fashion as the disk flux. We consider possible alternatives for the origin of the line. The connection between solar activity and 
the strength of Saturn's disk X-ray emission is investigated and compared with that of Jupiter. We also discuss the apparent lack of 
X-ray aurorae on Saturn; by comparing the planet's parameters relevant to aurora production with those of Jupiter we conclude that 
Saturnian X-ray aurorae are likely to have gone undetected because they are below the sensitivity threshold of current Earth-bound 
observatories. A similar comparison for Uranus and Neptune leads to the same disappointing conclusion, which is likely to hold true 
also with the planned next generation International X-ray Observatory. The next step in advancing this research can only be realised 
with in-situ X-ray observations at the planets. 

Key words. Planets and satellites: individual: Saturn 

1 . Introduction of Earth are solar wind driven (Clarke et al. l2005l l: Saturn's UV 

aurorae appear to be at variance with both these scenarios, rather 
In recent years X-ray observations have provided a novel way of than bd intermediate bet ween the two, as originally expected 
furthering the study of planets m our solar system and our un- (Bhardwaj and Gladstone [2000}. i n Saturn's case the UV au- 
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Early observations of Saturn with XMM-Newton and 
Chandra (Ness et al. 2004a, Ness et al. 2004b) unambiguously 
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Table 1. XMM-Newton and Chandra observations of Saturn 





Observation 


Observation 
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Geocentric 
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References 






(ks) 




(AU) 


(AU) 




XMM-Newton (1) 


2002 Oct. 01, 14:00 UT 


21 


18.8" 


9.0 


8.8 


Ness et al. 2004a. this paper 


Chandra (1) 


2003 Apr. 14, 18:00 UT 


66 


17.5" 


9.0 


9.5 


Nessetal.l2004bl 


Chandra (2) 


2004 Jan. 20, 05:30 UT 


37 


20.5" 


9.0 


8.1 


Bhardwai et al.|2005all2005bl 














this paper 


Chandra (3) 


2004 Jan. 26, 20:00 UT 


36 


20.4" 


9.0 


8.2 


Bhardwai et al.l2005all2005bl 


XMM-Newton (2) 


2005 Apr. 22, 10:00 UT 


84 


18.1" 


9.1 


9.2 


This paper 


XMM-Newton (3) 


2005 Oct. 28, 21:00 UT 


74 


18.2" 


9.1 


9.1 


This paper 



detected X-ray emission, found it to be concentrated in the equa- 
torial regions of the planet, and indicated that any X-ray auroral 
component, if present at all, is not easily detectable. More re- 
cent Chandra observations have provided us with a remarkable 
view of the X-ray morphology of the Saturnian system, and of 
its temporal behaviour in the X-rays: both the planet's disk and 
the rings have been separately detected, and strong variability, 
closely following the behaviour of the solar X-ray flux, has been 
discovered in the disk X-ray emission (Bhardwaj et al. 2005a, 
Bhardwaj et al. |2005bl >. This suggests that its origin may lie in 
the scattering of solar X-rays in the upper layers of Saturn's at- 
mosphere, as is the case for Jupiter's disk emission (Bhardwaj 
et al. |2005cl Bhardwaj et al. [20061 Cravens et al. 120061 Unlike 
Jupiter, though, no difference is found in the spectral character- 
istics of Saturn's disk and polar emissions, thus no evidence for 
X-rays of auroral origin is provided by the Chandra data. In fact, 
Bhardwaj et al. (2005a) remark that the X-ray emission from the 
polar regions appears to be anticorrelated with the strength of 
the FUV aurora, measured by the HST Imaging Spectrograph 
simultaneously with the Chandra observations. Moreover, only 
relatively little detail has been published so far about Saturn's 
X-ray spectral characteristics: an optically thin, coronal model 
is a good match for the disk emission (Ness et al. 12004a Ness 
et al. I2004bl Bhardwaj [2006 ), while the rings spectrum is prac- 
tically made up of a single O Ka emission line, which has been 
attributed largely to scattering of solar X-rays on the icy water 
ring material (Bhardwaj et al. 2005b). However, solar fluores- 
cence alone could not explain the ring X-ray brightness com- 
pletely. To follow up on these issues, we made more extensive 
XMM-Newton observations of Saturn in 2005: this was a particu- 
larly appropriate time to take a deeper look at the planet, while a 
campaign of HST observations was taking place, and the Cassini 
spacecraft (which entered orbit around Saturn in July 2004) was 
making coordinated in-situ measurements of a variety of physi- 
cal parameters in its environment. 

In this paper we report the results from our 2005 observa- 
tions, together with a re-analysis of the earlier (2002) data ex- 
tracted from the XMM-Newton archive (sees 2 and 3). We also 
model the Chandra X-ray spectra from the 2003 and 2004 obser- 
vations (which have already been described by Ness et al. 2004b, 
Bhardwaj et al. 2005a and Bhardwaj et al. 2005b - see sec. 3). 
We discuss our findings in sec. 4, also making comparisons with 
Jupiter and Earth, as well as Uranus and Neptune, and draw our 
conclusions in sec. 5. 

2. XMM-Newton observations of Saturn 

In Table 1 we list in time order all the X-ray observations of 
Saturn which have been carried out so far with XMM-Newton 



and Chandra, identifying them by sequential numbers; those 
with Chandra were all performed with the Advanced CCD 
Imaging Spectrometer (ACIS), so that spectral information was 
available. Table 1 includes the mid-epochs and durations of the 
observations, as well as the apparent size of the planet's disk, 
its distance from the Sun and the Earth at those times, and rele- 
vant references. During our 2005 XMM-Newton observations the 
CCD European Photon Imaging Cameras (EPIC-pn, Striider et 
al. l200T1 EPIC-MOS, Turner et al. l200Tb were operated in Full 
Frame mode with the thick filter at all times, to avoid any op- 
tical leak contamination (see Ness et al. 2004a for a discussion 
of the optical light leak affecting thin and medium filters). In 
2002 only one exposure (with the EPIC-pn camera) was carried 
out with the thick filter, so the analysis had to be restricted to it 
for this epoch. In order to make a direct and consistent compar- 
ison of all XMM-Newton data, the analysis of the 2005 obser- 
vations reported here is also limited to the EPIC-pn camera. We 
note, though, that the best fit of the combined EPIC-pn and MOS 
datasets for 2005 is consistent, within the parameters errors, with 
that obtained using the EPIC-pn camera on its own. During all 
XMM-Newton observations the Reflection Grating Spectrometer 
(RGS, den Herder et al. 2001 ) was run in Spectroscopy mode, 
but the source brightness is too low to enable high resolution 
spectroscopy to be carried out. The Optical Monitor telescope 
(Mason et al. [200Tb had its filter wheel kept in the BLOCKED 
position, because Saturn's optical brightness exceeds the safe 
limit for the instrument (thus no OM data were collected). 

At all 2002 and 2005 observing epochs Saturns's path on the 
sky was essentially in RA, and very slow (8"/hr max). Each one 
of the longer, 2005 observations was broken into two intervals, 
with a small pointing trim in between, to ensure the planet would 
be imaged on the same CCD chip at all times. The two intervals 
were combined during the data processing and the subsequent 
analysis. 

The two XMM-Newton 2005 observations, as well as the 
archival one from 2002, were processed and analysed in a uni- 
form fashion under SAS (Science Analysis Software) v. 7.1; the 
2002 observation (Ness et al. 2004a) had never been processed 
in this now standard way, because of the lack, at the time, of 
dedicated software to re-register the XMM-Newton events in the 
rest frame of a moving target. After this was done for the three 
Saturn datasets, images at the location of the planet were con- 
structed with the data from the EPIC cameras. Fig. 1 shows the 
EPIC-pn images of Saturn in the energy range 0.3 - 0.8 keV for 
the XMM-Newton observations 1 (top of the colour panels), 2 
(middle) and 3 (bottom; the planet is undetected on this occa- 
sion). Overposed to each image, at the nominal location of the 
planet's centre, is a white circle of 26" radius. At the top of Fig. 
1 is a diagram showing the appearence of Saturn and its rings in 



G. Branduardi-Raymont et al.: X-rays from Saturn: A study with XMM-Newton and Chandra over the years 2002-05 



3 



s 




5 



» r. , . i i .... i .... i .... i - 

5 55 32 5 55 31 5 55 30 5 55M 5 55 Z6 

Right Ascension <h m s) 




R 5 



Fig. 1. The black and white diagram in the top panel shows the 
appearence of Saturn and its rings (outer ellipse major axis = 
44") at the epoch of the XMM-Newton observations in Oct. 2002 
(the view in 2005 was very similar). In the three colour panels 
are XMM-Newton EPIC-pn images of Saturn in the 0.3 - 0.8 
keV band for Oct. 2002 (top), Apr.2005 (middle) and Oct. 2005 
(bottom). Overplotted in white are 26" radius circles centred on 
the planet, delimiting the area used for spectral extraction. The 
size of the diagram and of the images is 70" x 70", and the EPIC- 
pn image pixels are of 2.9" side (North is at the top and East to 



Oct. 2002 (the view, with the South Pole unobstructed, was very 
similar in Apr. and Oct. 2005). 

From visual comparison of the three panels in Fig. 1 it is 
clear that the planet was X-ray brightest in 2002, dimmer in Apr. 
2005 and weakest, in fact undetected, in Oct. 2005 (see sec. 3 
and Table 2 for X-ray upper limits at this epoch). 

Despite the relatively broad XMM-Newton Point Spread 
Function (PSF, of 15" Half Energy Width) compared with the 
planet's image size (~20" diameter), the morphology of the X- 
ray emission at both 2002 and Apr. 2005 epochs does not appear 
circularly symmetric, but is seen to stretch along the direction 
of the rings plane. In the following we will analyse the XMM- 
Newton data as encompassing both, disk and ring emissions. 

We also searched for evidence of X-ray emission from 
Saturn's satellite Titan in the XMM-Newton datasets, but none 
was found. Similarly no detection was obtained with the 
Chandra observations. 

3. Spectral fitting: XMM-Newton EPIC and Chandra 
ACIS data 

In principle Saturn occults the cosmic X-ray background origi- 
nating beyond it, so we could expect that background subtraction 
may not be required in the analysis of its spectra; however, the 
contribution of particle background (dominant at higher ener- 
gies), and any 'spill over' of the cosmic X-ray background onto 
the planet's disk due to the width of XMM-Newton PSF may be- 
come significant if we use a region to extract the X-ray events 
for the planet which is large compared with the size of its disk. 
This is our case, where we are using an extraction circle of 26" 
radius, in order to include as much as possible of the X-ray flux 
from the source (Saturn's disk diameter 20"). For an XMM- 
Newton point source (which, for this purpose, Saturn can be ap- 
proximated to), an extraction radius of 26" corresponds to an 
encircled energy of 80%, so all fluxes quoted in the following 
for the planet are corrected for this factor. Fig. 2 shows a com- 
parison of Saturn's 2002 and Apr. 2005 EPIC-pn spectra (black 
crosses) extracted from the regions delimited by white circles in 
Fig. 1 , with those (red crosses) from background regions (devoid 
of sources) of the same shape and size, from the same EPIC-pn 
CCD. The plots indicate that the emission from the planet lies 
above the background in the band 0.3 - 2 keV, and that it prac- 
tically disappears under it at higher energies. The XMM-Newton 
spectral fitting results presented below refer to the Oct. 2002 and 
Apr. 2005 observations (the planet was undetected in Oct. 2005). 

In the case of the Chandra data, unlike the analysis carried 
out by Bhardwaj et al. (2005a,b), we extracted the full planet 
data, i.e. the combined emission from the disk, the polar cap 
and the rings, for the purpose of comparing with XMM-Newton 
which cannot resolve these morphological details. As already 
pointed out by Bhardwaj et al. (2005a), the spectrum of Saturn's 
polar region is very similar to that of the planet's disk. Here we 
present results from the fits made to the Apr. 2003, the 20 and 
the 26 Jan. 2004 Chandra data. The 2003 observation was made 
up of two exposures: we show the best fits for the two datasets 
combined, since individually they give very similar results. The 
source was weaker by a factor of three on 26 Jan. 2004 than a 
week earlier (Bhardwaj et al. |2005a| l, when Saturn brightened 
in the X-rays in direct response to a powerful X-ray flare taking 
place on the solar hemisphere visible from the planet. 

Before fitting, the spectra were binned so as to include at 
least 15 counts per energy channel (10 counts for the low-flux 
Chandra observation of 26 Jan. 2004), to validate the applicabil- 
ity of x 1 testing for parameter estimation (Cash 1979). We used 
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Fig. 2. XMM-Newton EPIC-pn spectra of Saturn (black crosses, 
before any background subtraction) extracted from the regions 
marked by white circles in Fig. 1 ; in red are the spectra of the 
background, extracted from regions devoid of sources, of the 
same circular shape and size as for Saturn, and from the same 
EPIC-pn CCD (top: Oct. 2002; bottom: Apr. 2005). 



XSPEC v. 12.4.0 for the fitting and modeled the spectra with an 
optically thin mekal component (to represent the scattered solar 
X-ray contribution) and a Gaussian line at 0.53 keV to take into 
account the O Ka line emission from the rings. The mekal tem- 
perature and normalisation, and the energy and normalisation of 
the emission line were left free to vary in the fit. Elemental abun- 
dances for the coronal model were set at the solar level; the width 
of the emission line was kept fixed, and taken to be due only 
to instrumental broadening. This model combination provides 
good fits to the EPIC-pn spectra from the XMM-Newton obser- 
vations in 2002 and Apr. 2005, and also to the Chandra ACIS 
spectra for Saturn and its rings combined, from the 2003 and 
2004 observations. Fig. 3 shows the data (black crosses) and the 
best fits (red histograms) for all the XMM-Newton and Chandra 
datasets in time sequence (from top to bottom). 

Table 2 lists the total numbers of source counts in the spectra, 
the best fit parameters and their 90% confidence errors, as well as 
energy fluxes for the two spectral components separately, for all 
the spectra from the two observatories. The fluxes were obtained 
from the best fit models by setting to zero the normalisation of 
one or the other of the components and have been converted to 
emitted powers (also listed in Table 2) by multiplying by And 2 , 
where d is the geocentric distance of Saturn at each epoch. 

We have checked whether the presence of the oxygen emis- 
sion line in the spectra is statistically significant by setting its 
normalisation to and re-fitting the spectra. The resulting best 
fit^f 2 values are listed in the last column of Table 2 and demon- 
strate that ignoring the presence of the line worsens the fits sig- 
nificantly. 
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Fig. 3. XMM-Newton EPIC-pn and Chandra ACIS spectra of 
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Both Apr. 2003 and Apr. 2005 spectra in Fig. 3 show one bin 
deviating in excess of the best fit at ~0.75 keV (16.5 A); in both 
cases this feature can be accounted for by a narrow emission 
line at the instrumental resolution. Interestingly, Fe XVII lines 
at 15 and 17 A are characteristic of the solar coronal spectrum; 
however, the statistical quality of the data (and in particular the 
hard-to-constrain strength of the line) does not warrant exploring 
this further. 

The XMM-Newton (3) observation (Oct. 2005) shows no ev- 
idence of X-ray emission from the planet (see bottom panel of 
Fig. 1), and we extract a 3<x upper limit of 44 counts in the EPIC- 
pn camera over the full duration (74 ks) of the dataset. Adopting 
the spectral parameters that best fit the Apr. 2005 XMM-Newton 
(2) EPIC-pn observation, assuming the same relative strengths 
for the disk and line emissions, and normalising to the shorter 
exposure and the slightly smaller geocentric distance, we obtain 
an upper limit of 19 MW to the disk emitted power in Oct. 2005, 
down by more than a factor of 10 from the level observed three 
years earlier during the XMM-Newton (1) observation (Table 2). 

4. Discussion 

4.1. Separating Saturn's disk and rings X-ray emission 

Our spectral fitting of Saturn's XMM-Newton and Chandra data, 
gathered over the period Oct. 2002 - Oct. 2005, shows that the 
X-ray emission is well described by a combination of two com- 
ponents, an optically thin coronal model (mekal) and an O Ka 
emission line; in accordance with the earlier results of Ness et al. 
(2004a,b) and Bhardwaj et al. d2005al) . we interpret the mekal 
component in Saturn's X-ray spectrum as representing its disk 
emission, which varies over the years 2002 - 2005 according to 
the diminishing strength of the solar X-ray output. The average 
temperature of the coronal emission is 0.5 keV, which is similar 
to that measured (0.4 - 0.5 keV) on three occasions in 2003 for 
Jupiter's disk X-rays (Branduardi-Raymont et al. [2007a); this 
strengthens the idea that in both cases the emission originates 
from the scattering of solar X-rays in the planets' upper atmo- 
spheres (cf. Bhardwaj et al. 2005a,c, Bhardwaj 2006, Cravens et 
al. l2006l 

The presence of an additional spectral component, in the 
form of an oxygen emission line, had already been suggested by 
the early fits of the 2002 and 2003 XMM-Newton and Chandra 
data by Ness et al. (2004a,b); this line was later associated 
with the planet's rings by Bhardwaj et al. (2005b) using the 
Chandra data, by spatially selecting only events coincident with 
the Saturnian rings. We have taken the different approach of at- 
tempting to separate spectrally the two emission components 
through a uniform analysis of the spectra of the whole planet 
+ rings system from both observatories. This solves the problem 
of the spatial resolution of XMM-Newton being too limited to al- 
low us to separate spatially the planet's disk from the rings; then 
we can study the variability of both emission components over 
the solar cycle by examining XMM-Newton and Chandra spec- 
tral fitting results, under the assumption that all the line emission 
originates in the rings. 

As a consistency check, we have compared the parameters of 
our best fits and the values of emitted power derived from them 
(Table 2) with those previously published for the 2002 XMM- 
Newton observation and the 2003 and 2004 Chandra datasets. In 
the first case, Ness et al. ( j2004a| ) use the apec model in XSPEC 
to describe the scattered solar spectrum, and combine it with a 
narrow oxygen line: their best fit temperature of 0.33 keV is con- 
sistent with our value of 0.28 keV (for a mekal model) within 
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Fig. 4. Saturn's power emitted in the oxygen line plotted against 
the 0.2 - 2.0 keV power from the planet's disk, for all XMM- 
Newton and Chandra observations. 

the errors; they derive apec and line fluxes of 1.36 and 0.22 x 
10~ 14 erg cirT 2 s _1 respectively, which are 30% and 10% higher 
than our mekal and line fluxes respectively, but again consis- 
tent within the spectral parameters errors. For the Chandra 2003 
observation, Ness et al. (2004b) fit only the emission from the 
circular region coinciding with the planet's disk (30% of which 
is occulted by the rings), deriving a best fit mekal temperature of 
0.39 keV, with a 0. 1 - 2 keV flux of 0.6 X 10~ 14 erg cirT 2 s~\ and 
a line flux of 0.13 x 10~ 14 erg cirT 2 s _1 . In this case our results 
are somewhat different (see Table 2) because we include the full 
rings emission in the fit (thus our oxygen line is stronger). We 
can make a more direct comparison with the results obtained by 
Bhardwaj et al. (2005b) on the rings line emission observed by 
Chandra. These authors quote a line flux of 0.5 x 10 erg cm -2 
s for the 20 Jan. 2004 observation which is in very good agree- 
ment with our value (0.48 x 10~ 14 erg cirT 2 s~'). Note that their 
flux is derived by fitting the rings spectrum on its own in XSPEC, 
while here we obtain the rings flux by extracting and modelling 
the combined disk and rings spectrum (in order to compare with 
the un-separated spectra from XMM-Newton). The consistency 
of the two results gives confidence in our approach. However, 
the line powers which we derive from the 2003 and 26 Jan. 2004 
Chandra data are factors of ~ 2 and 1.7 above what Bhardwaj et 
al. (2005b) quote respectively. A possible explanation of the dis- 
crepancy may be the presence of oxygen line emission from the 
rings area occulting the disk, which was excluded by Bhardwaj 
et al. but is included in our analysis. In fact, Ness et al. (j2004b) 
find that the inclusion of a small amount of line emission does 
improve their fit of the Chandra 2003 full disk region (including 
the part covered by the rings). If this is the case, the additional 
line contribution from this rings area may become more relevant 
when the disk flux is lower, as indeed can be seen in the Chandra 
observations (1) and (3) (Fig. 3). 

Turning to the power emitted by the disk, inspection of 
Bhardwaj et al. (2005a) Fig. 4 indicates 0.2 - 2.0 keV powers 
of -380, 170, 250 and 90 MW in 2002, 2003 and on 20 and 
26 Jan. 2004 respectively; these are a factor between 1.5 and 2.4 
larger than those we derive and list in Table 2. The different anal- 
ysis approach may be responsible for this discrepancy. However, 
since we are interested in studying trends over time, our uniform 
approach to the whole dataset ought to give more comparable 
results. 
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0.15+0.06 


82« 8 
74+31 

36+15 


12.2/10 

0.22/1 

17.4/16 


16.7/12 
5.2/3 
33.3/18 


3.23 
0.16 
0.22 


2005 


XMM-Newton (3) 


<44" 


0.33 


< 0.3'' 


< 0.08* 


< 19'' 


0.52 


< 0.51'' 


< 0.05'' 


< 12'' 






0.004 



" mekal temperature in keV 

b mekal normalisation at 1 keV in units of 10~ 6 ph cm" 2 s _I keV -1 

c Total 0.2 - 2.0 keV energy flux in the mekal component in units of 10 -14 erg cirT 2 s , the errors on the fluxes were obtained by varying the 

mekal temperatures and normalisations between their 90% confidence limits 
d Powers and errors calculated by multiplying the values in the previous column by 47rcP, where d is the geocentric distance of Saturn at the 

given epoch. 

e Gaussian line normalisation in units of 10~ 6 ph cirr 2 s -1 
' Total energy flux in the line in units of 10~ 14 erg cirr 2 s -1 
s x 2 value and degrees of freedom 
'' 3cr upper limit 

' Solar X-ray flux (1.5 - 12.4 keV) in units of 10~ 6 Watt nT 2 , corrected for solar rotation (see text) 



Table 3. Planetary parameters relevant to auroral studies 



Planet 


Average solar 


Rotation 


Equatorial 


Surface 


Magnetic 


FUV aurora 


X-ray aurora 


Disk X-ray 


Magnetospheric 




distance 


period 


radius 


magnetic field 


moment 


average input 


total average 


average 


ions max. 




(AU)" 


(hrs) n 


(Earth = Vf* 


(G) fl 


(Earth = If 


power (GW) f 


power (MW) rf 


power (MW) 


density (cirr 3 )" 


Earth 


1.0 


23.9 


1 


0.31 


1 


1 - 100 


10-30 


40 rf 


1 - 4000 


Jupiter 


5.2 


9.9 


11.2 


4.28 


20,000 


few x 10 4 


400 - 1000 


500 - 2000 rf 


> 3000 


Saturn 


9.6 


10.7 


9.4 


0.22 


600 


10 2 - 10 3 


7 


< 230 e 


100 ? 


Uranus 


19.2 


17.2 


4.0 


0.23 


50 


<40f 


? 


? 


3 


Neptune 


30.1 


16.1 


3.9 


0.14 


25 


<0A f 


? 


? 


2 



° From The New Solar System, Ed.s J. K. Beatty, C. Collins Petersen & A. Chaikin, 1999, Cambridge University Press 
* For the outer planets, radius at 1 bar atmospheric pressure 
c Clarke et al. ( f2005t 
d Bhardwaj et al. i2007l 
e This paper 

f Bhardwaj and Gladstone ( [20001 
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Then we can use the whole sequence of XMM-Newton and 
Chandra observations to study how the fluxes in the two X-ray 
spectral components, the emission from the planet's disk and the 
oxygen line (assumed to be confined to the rings), varied rela- 
tive to each other over a period of three years. From Fig. 3 it 
is already clear that the line is more prominent when the disk 
flux is lower (Apr. 2003, 26 Jan. 2004 and Apr. 2005) and is less 
visible when the disk is brighter (Oct. 2002 and 20 Jan. 2004). 
This suggests that the two components vary over time in different 
ways: Fig. 4 shows the power emitted in the oxygen line plotted 
against that from the disk (from Table 2). The trend of the disk 
emission to decrease with time (mirroring the decay of the solar 
X-ray output as the solar cycle progresses towards a minimum), 
originally identified by Bhardwaj et al. (2005a), is apparent from 
the temporal sequence of the data in Fig. 4. However, despite the 
large errors, it is clear that there is no similar trend in the oxygen 
line; although both disk and line powers change by factors of ~ 6 
and 4 respectively over the three years of observations, their vari- 
ability appears to be un-correlated: in particular, the line strength 
is consistent with being the same in 2002 and Apr. 2005, while 
the disk has dimmed by a factor of 3.5. In addition, as reported 
by Bhardwaj et al. (2005b), the response of the disk to the solar 
flare on 20 Jan. 2004 is not clearly matched by a similar increase 
in the rings line power. Both disk and rings have disappeared by 
Oct. 2005, which may argue for some correlation between the 
two; however, it is not clear that the line, if present but weak, 
would produce enough counts to be detectable on its own. 

4.2. The solar connection 

Fig. 5 shows Saturn's disk 0.2-2.0 keV X-ray power from Table 
2 (open squares) plotted against the solar 1 - 8 A (1.5 - 12.4 
keV) flux (also listed in Table 2) measured by the GOES-10 
satellite at each observing epoch and averaged over the dura- 
tion of the observations. These fluxes have been corrected for 
the Sun's rotation except in the cases (Jan. 2004) when the Earth 
and Saturn were practically in opposition. We have neglected the 
travel time differences Sun-Earth and Sun-Saturn-Earth because 
these are very small with respect to the correction for solar rota- 
tion (and cancel out in opposition). As expected, there appears to 
be correlation between Saturn's disk X-ray power and the solar 
X-ray flux. For comparison, the same parameters are also plotted 
for Jupiter's disk (filled circles), which also has been shown to 
act as a mirror of solar X-rays (Maurellis et al. 2000 Bhardwaj 
et al. |2005cl Cravens et al. 2006). XMM-Newton and Chandra 
observations in Fig. 5 are identified by an X and a C, respec- 
tively. 

For completeness here is a brief summary of Jupiter's mea- 
surements used in Fig. 5. Waite et al. (1997) were the first to 
report an X-ray measurement for the low latitude regions of the 
planet from their ROSAT High Resolution Imager observation in 
July 1994. Chandra observed Jupiter with the High Resolution 
Camera (HRC) in Dec. 2000, with both the HRC and ACIS in 
Feb. 2003 and with ACIS in Feb. - March 2007; this last time 
the observations were simultaneous with XMM-Newton, which 
had looked at the planet also in April and Nov. 2003. Bhardwaj 
et al. ( 12006b report 0.5 - 1.5 keV fluxes for the two 2003 ACIS 
observations, which we have converted to 0.2 - 2.0 keV using 
their mekal best fit parameters. They also mention that the HRC 
flux was about 50% of that in Dec. 2000, so we have attributed 
twice the mean power measured by ACIS in 2003 to the Dec. 
2000 epoch (when Gladstone et al. [20021 only had HRC data, 
without spectral information with which to characterise the spee- 
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Fig. 5. Saturn's (open squares) and Jupiter's (filled circles) disk 
0.2 - 2.0 keV X-ray power plotted against the solar 1 - 8 A (1 .5 
- 12.4 keV) X-ray flux measured by the GOES-10 satellite, for 
all XMM-Newton (X) and Chandra (C) observations. The 1994 
Jupiter data came from ROSAT. The continuous lines are linear 
square fits to each dataset. 



trum properly). Emitted powers for the XMM-Newton observa- 
tions of Jupiter in 2003 are reported in Branduardi-Raymont et 
al. (2007a), but were computed using a factor of 2n between 
flux and power, so we doubled them to be comparable with the 
rest of the measurements. The 2007 XMM-Newton observations, 
which will be the subject of a future publication (Branduardi- 
Raymont et al., in preparation), took place at the time of the New 
Horizon Jupiter fly-by, and were contemporaneous with some by 
Chandra; at this epoch, close to solar minimum, the disk X-ray 
power had dimmed considerably (~120 MW). The errors asso- 
ciated with the 2003 XMM-Newton mekal model normalisations 
are of the order of 10% (Branduardi-Raymont et al. 20073, so 
we attributed the same fractional error to the corresponding emit- 
ted powers, and to the Chandra measurements for which no er- 
rors are quoted. All these errors are contained within the size of 
the plotting symbols. Larger errors (+100%, -50%) apply to the 
1994 ROSAT power (Waite et al \TWj\ . The linear correlation 
coefficients are 0.60 and 0.31 for the Saturn and Jupiter datasets, 
respectively. The correlation in Jupiter's case becomes stronger 
if we exclude the 1994 and 2000 measurements, for which we do 
not have spectral information (the coefficient increases to 0.54). 
Plotted in Fig. 5 are also lines derived from a linear square fit to 
each dataset, which reinforce the view of the flux trends. The 
correlation between Saturn's and Jupiter's disk X-ray powers 
and the GOES 0.5 - 3 A (4 - 25 keV) flux is much weaker (co- 
efficients of 0.39 and 0.18 respectively). This is not surprising, 
since the band covers higher energies than those at which the 
~0.5 keV solar spectrum, scattered by the planet's atmosphere, 
peaks. 

4.3. The origin of the rings X-ray emission 

Our analysis of the complete set of Saturn's X-ray spectra over 
the three year period from Oct. 2002 to Oct. 2005 suggests that 
the strength of the oxygen emission line from its rings does not 
vary in the same way as the disk emission does, following so- 
lar activity. This would be at variance with the interpretation of 
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the origin of the rings emission as fluorescent scattering of so- 
lar X-rays from the icy rings material (Bhardwaj et al. [2005b). 
However, only about one third of the ring X-ray emission could 
be explained in this way by Bhardwaj et al. (2005b). While the 
fluorescence interpretation is attractive because is underpinned 
by the same mechanism which leads to the disk X-ray emission, 
and our spectral separation technique cannot prove that the line 
emission is truly all confined to the rings, nevertheless it is in- 
teresting to explore a possible alternative to the solar scattering 
interpretation. 

Chandra imaging in both 2003 and 2004 appears to indicate 
that the X-ray emission concentrates on the East ansa (morning 
side) of the rings, where also spokes occur; indeed Bhardwaj 
et al. (2005b) suggest that the higher X-ray brightness could be 
due to enhanced fluorescence on the ice-dust particles that are 
thought to make up the spokes, and to be created by meteoritic 
impacts, which are more likely in the early morning hours. On 
the other hand, Cassini measurements have shown the spokes to 
be related to electrons injection (Coates et al. 2005). Recently it 
has also been suggested that lightning-induced electron beams 
(following cosmic rays impacts) may produce the spokes (Jones 
et al. 120061 ). Such electron beams may also energise the icy par- 
ticles in the rings, which could then fluoresce and cause the oxy- 
gen line emission. Thus, if the X-ray line is related in this way 
to thunderstorms that are characteristically 'seasonal', since sea- 
sons last a long time on Saturn (~30 year sidereal period), one 
could expect the line strength to be more or less constant over 
the three year period for which we have observations. 

On the other hand, recent modeling work by Tseng et al. 
(2009) suggests that the oxygen atmosphere of the rings is 
strongly dependent on Saturn's position in its orbit, and a much 
larger density is expected in the Southern rather than in the 
Northern hemisphere. Thus it should not come as a surprise that 
the oxygen X-ray emission from the rings may vary over dif- 
ferent timescales from those of the disk emission, which is con- 
trolled by solar activity. Future observations may also be able 
to show whether the brightness of the oxygen X-ray emission 
changes after the equinox (11 August 2009). 

4.4. Are there X-ray aurorae on Saturn? and on Uranus and 
Neptune? 

Unlike for Jupiter (Branduardi-Raymont et al. 2007b), there is 
no evidence of auroral X-ray emission in the Saturn data ob- 
tained so far. The X-ray spectrum of the Jovian aurorae is made 
up of two components: the softer one is generally assumed to be 
produced by energetic heavy ions precipitating in the planet's at- 
mosphere and radiating line emission following charge exchange 
with atmospheric particles; the other, becoming dominant above 
2 keV, is due to energetic electrons also precipitating through 
the atmosphere and radiating bremsstrahlung. Recently, a mor- 
phological study of the spatial distribution of this high energy 
emission in Jupiter's Northern polar region has shown it to be 
co-located with the main FUV auroral oval and bright FUV re- 
gions (Branduardi-Raymont et al. 2008). This, and the relative 
emission powers in the FUV and hard X-rays, support the view 
that on Jupiter these two emissions are produced by the same 
population of energetic electrons precipitating through the at- 
mosphere. Given the bright and variable FUV aurorae observed 
on Saturn, one would expect some correlation with the X-ray 
emission to be present. 

On the basis of model calculations by Singhal et al. d 19921 ) 
and Waite et al. d 19921 ), Branduardi-Raymont et al. ( 120081) de- 
duce that the power in the FUV emission is expected to be some 



10 5 to 10 6 times that in hard X-rays if the same electrons are 
producing both and conclude that the observed emissions from 
Jupiter's aurorae (i.e. 45 MW in hard X-rays and 340 GW in 
FUV) are broadly consistent with these ratios. We can make sim- 
ilar considerations for Saturn, which also possesses a dense at- 
mosphere, like Jupiter's: restricting the analysis to the more sen- 
sitive XMM-Newton observations, we extract 3cr upper limits in 
the band 2.0 - 8.0 keV of 16, 35 and 37 counts for datasets (1), 
(2) and (3), respectively. Adopting a bremsstrahlung temperature 
of 10 keV (lower than the 90 keV best fitting the Jupiter spec- 
tra of Branduardi-Raymont et al. 2007b, since the precipitating 
electrons are generally of lower energy on Saturn than on Jupiter, 
Gustin et al. [20091 Gerard et al |20"U9l ), and using the EPIC-pn 
response matrices to determine the conversion factor from coun- 
trate to energy flux, we find corresponding 3<x upper limits to 
the emitted 2.0 - 8.0 keV power of 147, 94 and 107 MW. The 
XMM-Newton (3) observation (28 Oct. 2005, 06:34 to 29 Oct. 
2005, 11:32 UT) was contemporaneus with three exposures of 
the Saturnian aurora taken with the HST Advanced Camera for 
Surveys, during a coordinated campaign with in-situ measure- 
ments by the Cassini spacecraft. This was a time of quiet solar 
wind conditions, and yet the UV aurora was seen to be very vari- 
able; the measured UV power at the time of the XMM-Newton 
observation varied between 3.3 and 8.0 GW (Gerard et al. 2006). 
If we apply to Saturn the ratios of 10 5 to 10 6 between hard X- 
ray and FUV powers calculated for Jupiter, it is clear that any 
electron bremsstrahlung emission would have been some 3 to 4 
orders of magnitude fainter than our current XMM-Newton upper 
limits. 

To make a further comparison of the most magnetised plan- 
ets in the solar system, in Table 3 we summarise parameters that 
can be expected to be relevant to the auroral characteristics for 
all the outer planets, as well as Earth. With respect to Jupiter, 
Saturn is down by a factor of ~ 20 and 30 in surface magnetic 
field and magnetic moment respectively. A similar ratio applies 
to the maximum average FUV aurora input power. Moreover, 
other parameters are likely to matter in the aurora production, 
such as the magnetospheric particle density: this is significantly 
higher, by at least another factor of 30, around Jupiter (being 
fed by Io's volcanoes) than Saturn (where it is sustained by the 
rings and the icy moons). Thus, if Saturn's total auroral X-ray 
power (including both soft X-rays of ionic origin and the higher 
energy electron bremsstrahlung) is lower than that of Jupiter by 
the combination of the magnetic field and particle density factors 
(around three orders of magnitude), it may measure only a few 
MW; this is equivalent to an X-ray flux at Earth of a few x 10" 17 
erg cm" 2 s" 1 , so it is no surprise that it is going undetected with 
current instrumentation. In fact, Saturn's aurora may be at the 
limit of detectability even with the next generation most pow- 
erful International X-ray Observatory (IXO), being planned by 
ESA, NASA and JAXA, which is expected to have a sensitivity 
limit of ~ 10" 17 erg cm" 2 s" 1 . 

One may argue that our scaling approach is rather simplis- 
tic, but there are other examples in the planetary context where 
it seems to be applicable. Radio events reported by Louarn et al. 
(2007) at both Jupiter and Saturn suggest that plasma transport 
in the magnetosphere and the associated release of rotational en- 
ergy may operate in a similar way at the two planets, but scale 
according to different timescales for internal processes. This is 
also in line with the conclusions of Stallard et al. (2008) that 
the differences in auroral emissions between the two giant outer 
planets are determined by scaling differences in their physical 
conditions and environments, while the processes leading to au- 
rora formation are essentially the same (e.g. breakdown of co- 
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rotation). There are terrestrial analogues to the scaling idea as 
well: solar wind disturbances at Earth are larger in amplitude and 
faster (order of minutes) than at Jupiter and Saturn (timescales 
of hours) because of the smaller magnetosphere and the lack 
of steady supply of plasma (Clarke et al. 2009 and references 
therein). 

Assuming that the same considerations are valid for Uranus 
and Neptune, and scaling on their dipole magnetic moments, 
which are down by a factor of 400 and 800 from Jupiter re- 
spectively, we would expect FUV auroral strengths of the or- 
der of 100 to 50 GW in the two cases. However, the magne- 
tospheric particle density at these planets is some 10 3 times 
lower than at Jupiter (Table 3), suggesting that we should expect 
much lower FUV emitted powers than this initial extrapolation 
implies. The observed FUV powers, reported by Bhardwaj and 
Gladstone ( 120001 ), are at most 40 and 0.1 GW: the latter value is 
of the order which we would expect for the low plasma content 
of Neptune's magnetosphere, but Uranus does not seem to scale 
down in the same way; perhaps the very large tilt (59°) between 
its magnetic dipole and rotation axis and the 30% offset of the 
dipole from the planet's centre, causing an order of magnitude 
difference in surface magnetic field strength between day- and 
night-side, have also the effect of enhancing the power output in 
the aurora. In any case, taking these estimates together with the 
vaste distances from Earth of these two planets excludes the pos- 
sibility of detecting their X-ray auroral emissions with current 
instruments. Indeed Bhardwaj et al. (2007) report an unsuccess- 
ful attempt by Chandra to detect X-rays from Uranus in August 
2002. These authors also suggest that Uranus and Neptune, as 
well as Saturn's moon Titan, are likely to shine in X-rays, both 
by particle precipitation in the magnetosphere, and by scattering 
of solar X-rays. Looking to the future, let us consider the most 
optimistic case for Uranus, at about 4 times the distance from 
Earth than Jupiter, with an FUV output ~ 1000 times lower than 
Jupiter: ignoring the much smaller magnetospheric particle den- 
sity, and scaling from the observed Jovian 0.2 - 2 keV auroral 
emission (few times 10~ 14 ergcrrT 2 s _1 , Branduardi-Raymont et 
al. 12007b) implies a flux at Earth from Uranus of ~ 10~ 18 erg 
cm -2 s _1 , far too low for current X-ray observatories, and about 
an order of magnitude below IXO's expected sensitivity limit. 

5. Conclusions 

1. Statement 1. We report the results of all the X-ray obser- 
vations of Saturn carried out to date. We present a detailed 
spectral analysis of the available data which reinforces the 
earlier conclusion (Bhardwaj et al. |20Q5al > that the planet's 
disk acts as a reflector of solar X-rays, much in the same way 
as observed for Jupiter, Earth, Venus and Mars (cf. Bhardwaj 
etal. 12007b . 

2. Statement 2. Our results also suggest that the strength of the 
X-ray emission from the Saturnian rings, in the form of a flu- 
orescence oxygen line, does not systematically decrease, like 
that from the disk, while the solar activity cycle evolves to- 
wards its minimum. As a possible alternative to scattering of 
solar X-rays in the rings icy material, we suggest that the ori- 
gin of the line may be related to the formation of the spokes, 
which appear on the East ansa of the rings as do the bright 
X-ray spots seen by Chandra, by lightning-induced electron 
beams. 

3. Statement 3. Unlike Jupiter and Earth, we do not find ev- 
idence for X-ray aurorae on Saturn. This can be explained 
by the limited sensitivity of the observations to date, if the 
auroral X-ray and FUV emissions on Saturn have relative 



strengths similar to those on Jupiter. Other factors, such as 
magnetic field strength, magnetospheric particle densities 
and internal plasma sources may also contribute. 
4. Statement 4. After comparison with Jupiter and Earth, we 
conclude that Saturn is unique in its X-ray characteristics as 
it is in the FUV (Clarke et al. 2005). We still have much 
to learn about the precise mechanisms operating in its auro- 
rae, and in its rings, leading (or not, as it may be) to X-ray 
production. The effectiveness of further observations with 
Chandra and XMM-Newton will increase as we move out 
of the current (2009) minimum in the solar activity cycle, 
providing invaluable additional information on the evolution 
of both, disk and rings emissions, and of their response to 
the solar X-ray flux and the solar wind. However, a much 
more efficient way to pursue and enhance these studies, and 
relate the electromagnetic output of the planet to particle and 
magnetic field measurements, would involve flying an X-ray 
spectrometer to operate in-situ, onboard a spacecraft mission 
to the Saturnian system. We look forward to the possible re- 
alisation of this project in association with the missions to 
the outer planets currently being studied by both ESA and 
NASA. Such a mission to Saturn is most likely to involve a 
close encounter with Titan, and an X-ray spectrometer would 
return unique data on this enigmatic satellite, and its interac- 
tions with Saturn and the solar wind as well, especially at 
times when it is outside Saturn's magnetosphere. In-situ X- 
ray observations are also the only feasible way to search for 
auroral X-ray emission from Uranus and Neptune. 
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